arXiv: 1509.01803vl [cond-mat.quant-gas] 6 Sep 2015 


Traveling Majorana solitons in a one-dimensional spin-orbit coupled Fermi superfluid 

Peng Zou,^ Joachim Brand,^ Xia-Ji Liu,^ and Hui Hu^’0 

^ Centre for Quantum and Optieal Seienee, Swinburne University of Teehnology, Melbourne 3122, Australia 
^New Zealand Institute for Advaneed Study, Centre of Theoretieal Chemistry and Physies, 
and Dodd-Walls Centre for Photonie and Quantum Teehnologies, Massey University, Auekland, New Zealand 

(Dated: September 8 , 2015) 

We investigate traveling solitons of a one-dimensional spin-orbit coupled Fermi superfluid in both 
topologically trivial and non-trivial regimes by solving the static and time-dependent Bogoliubov- 
de Gennes equations. We find a critical velocity Vh for traveling solitons that is much smaller 
than the value predicted using the Landau criterion due to the presence of spin-orbit coupling, 
which strongly upshifts the energy level of the soliton-induced Andreev bound states towards the 
quasi-particle scattering continuum. Above Vh, our time-dependent simulations in harmonic traps 
indicate that traveling solitons decay by radiating sound waves. In the topological phase, we predict 
the existence of peculiar Majorana solitons, which host two Majorana fermions and feature a phase 
jump of TT across the soliton, irrespective of the velocity of travel. These unusual properties of 
Majorana solitons may open an alternative way to manipulate Majorana fermions for fault-tolerant 
topological quantum computations. 

PACS numbers: 03.75.Lm, 67.85.Lm, 67.85.De 


Solitons or localized waves that arise from the inter¬ 
play between the dispersion and nonlinearity of underly¬ 
ing systems are fascinatin^henomena occurring in many 
different fields of physics [l|. Over the past two decades, 
a major research emphasis has focused on solitons in 
atomic Bose-Einstein condensates (BECs) [H. The fam¬ 
ily of BEG solitons consists of many interesting members, 
from bright solitons in attractive BECs @ and gap soli¬ 
tons in optical lattices Q, to dark solitons in repulsively 
interacting BECs [3|-0 , which are created experimentally 
by imprinting a sharp and characteristic phase jump into 
the BEC. Remarkably, dark solitons may also be created 
in strongly interacting Eermi gases |8l4l6| at the crossover 
from BECs to Bardeen-Cooper-Schrieffer (BCS) super- 
fluids 0 , where phase kinks are encoded in the pairing 
order parameter. Their recent experimental observation 
may offer valuable insights into the nature of fermionic 
superfluidity in the strongly correlated regime 13, 18| . 


In this Letter, we consider traveling fermionic solitons 
in a different setup - one-dimensional (ID) Eermi super- 
fluids with spin-orbit coupling (see Eig. [T]) - and predict 
the existence of an exotic member of the soliton family 
when the superfluid becomes topologically non-trivial. It 
is referred to as Majorana soliton, owing to its ability 
to host two Majorana fermions that obey non-Abelian 
statistics at the soliton core 00. Majorana solitons 
are universal and remarkably robust, in the sense that 
their properties are not affected by a finite velocity of 
travel. In particular, the phase jump across a Majorana 
soliton is exactly pinned to tt and the density profile is 
unchanged. In other words, Majorana solitons are not 
greyed by a finite velocity. This unique stability renders 
Majorana solitons an ideal platform to manipulate Ma¬ 
jorana fermions for practical a ppl ications such as topo¬ 
logical quantum computations [21| . 
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Figure I: (color online). Upper panel: Sketch of the proposed 
experimental configuration. Many ID tubes of spin-orbit cou¬ 
pled Fermi gases are formed using a 2D optical lattice and 
two counter-propagating Raman lasers. Lower panel: The 
magnitude |A(^)| and phase 0(^) of the soliton order param¬ 
eter in the non-topological phase with parameters 7 = 3.41, 
h = 0.52Ef and Mf/Ef = 1.71. The solid, dashed and dot¬ 
ted lines correspond to soliton velocities Vs — 0 , 0A5vf and 
0.3 uf, respectively. 


Our investigation is motivated by the recent realiza¬ 


tions of spin-orbit coupling in atomic Fermi gases [22 


and the promising perspective of creating an atomic topo¬ 
logical superfluid 


24 l25j . Traveling Majorana solitons 
with fixed tt phase step, if experimentally observed to 
oscillate inside a Fermi cloud, would be a smoking-gun 
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proof of the existence of long-sought topological superflu¬ 
ids. We note that stationary dark solitons with Majorana 
fermions in a spin-orbit coupled Fermi gas were recently 
predicted [26|, [27|. However, the crucial issue raised in 
any practical manipulations, i.e., the fate of these soli- 
tons at a finite velocity of motion, was not addressed. 

Our results also suggest that a critical velocity for 
the stability of traveling solitons is greatly affected by 
spin-orbit coupling. This is because one of the mid-gap 
energy levels, the soliton-induced Andreev bound states 
(ABSs), is strongly up-shifted towards the bulk quasi¬ 
particle scattering continuum by spin-orbit coupling. As 
a result, the Andreev bound states are lost at a much 
smaller velocity Vh than the Landau critical velocity, 
which leads to a decay channel. At a velocity above Vh, 
we find that traveling solitons gradually decay via radi¬ 
ating sound waves. 

Model. We start by describing a possible experimental 
configuration, as sketched in the upper panel of Fig.[TJ A 
bundle of parallel, identical ID spin-1/2 Fermi gases 
can be formed by adding a tight 2D optical lattice in the 
transverse y — z plane [28|, and the spin-orbit coupling 
with equal Rashba and Dresselhaus weight can be real¬ 
ized by adapting the so-called NIST scheme using two 
counter-propagating Raman laser beams [22|. The re¬ 
sulting ID spin-orbit coupled Fermi gas in a single tube 
is modeled by the Hamiltonian H = f dx [Hol-Lint], 
where 


i’1- (x) " 

( 1 ) 

is the spin-orbit coupled single-particle part and 

'Hint = flioV’t (x) d (x) -Ipl (x) -t/jf (x) (2) 

with giD < 0 is the interaction Hamiltonian describing 
the attractive contact interaction between the two spin 
states (cr Here, 'ipl is the fermionic field opera¬ 

tor that creates an atom with mass m in the spin state 
cr. The term XkxCJy — hcTz with the momentum operator 
kx = —idjdx and Pauli matrices ay and az is induced by 
the Raman process, describing a synthetic spin-orbit cou¬ 
pling with strength X = fi^kRjm and an effective Zeeman 
field h = fli?/2, where kR and Qr are the momentum 
and Rabi frequency of the Raman beams [Hj, respec¬ 
tively. The term Hs = —k?d^/(2rri) + Vt{x) — y with the 
chemical potential fi describes the motion of atoms in a 
harmonic trapping potential Vt{x) = vraJ^x^ j2. 

We solve the model Hamiltonian for stationary and 
traveling solitons within the mean-field approxima¬ 
tion. This amounts to finding solutions with phase- 
twisted order parameter in the static and time-dependent 
Bogoliubov-de Gennes (BdG) equations, HBdG^yix) = 
Ey^r]{x) and HBdG^yixH) = ih{d/dt)^r]{xH), re¬ 
spectively. Here, for convenience we have used the 
Nambu spinor representation and have introduced = 


Ho = 


d (x), d (x) (Hs + \kx<Xy - ha^ 


wavc-functiou and en¬ 
ergy of Bogoliubov quasiparticles. The BdG Hamiltonian 
reads 


'^BdG = 


Hs — h —Xd/dx 
Xdjdx Hs + h 
0 A* 

-A* 0 


0 -A 

A 0 

—HsEh Xdjdx 

—Xdjdx —Hs — h 


( 3 ) 


and the BdG equations, either static or time-dependent, 
should be self-consistently solved with the gap equa- 
tion A = 

and the number equation n = —(1/2) '^y[\uari\‘^f {Eri) + 
\van\‘^f{-Er,)], where f{E) = 1/(1 + is the 

Fermi-Dirac distribution function and the summation is 
performed for the energy level (labeled by rj) up to a 
high-energy cut-off Ec^ i.e., \Ey \ < Ec. 

To obtain a moving soliton in a trapped gas, we first 
find a stationary dark soliton at xo away from the trap 
center [l^. By evolving such an initial state in time, the 
soliton is accelerated by the trap potential and caused 
to oscillate inside the Fermi cloud. The same procedure 
has previously been used to understand the dynamics of 
dark solitons in a BEC-BGS Fermi superfluid and 
could also be employed in experiment. 

We also search for traveling soliton solutions on 
a homogeneous (untrapped) background that satisfy 
A(x,t) = A{x — Vgt) = A(^), by solving the BdG eq na- 
tions in the co-moving frame with the velocity Vs |llj : 




Ey-ihvs-^ 




( 4 ) 


Here, HBdGiO obtained by replacing dx with d^ and 
A{xH) with A(^) in Eq. (|3j). In other words, we seek 
traveling solitons in a homogeneous gas that are station¬ 
ary in the frame of the soliton. This technique provides 
more insights into the soliton properties and enables us to 
isolate effects caused by the trappin g po tential when we 
analyze time-dependent simulations [l2|. Eor the calcu- 
ations in a box with length L we impose a modified peri¬ 
odic boundary condition, A(^ + L/2) = A(^ — Lj2)e'^^^, 
to explicitly take into account a phase jump Scf) across 
the soliton [l^ . In addition, we implement a generalized 
secant (Broyden’s) approach to make sure that the self- 
consistent iteration procedure will converge to a stable 
solution mi Hi. 

In numerical calculations, we use a dimensionless inter¬ 
action parameter to characterize the interaction strength, 
7 = —nng\Dj{jj^n\ which is basically the ratio between 
the interaction and kinetic energy at the density n. We 
choose the Eermi vector and energy, kp — 7rn/2 and 
Ep = h?k\j(2nnb)^ as the units of wave-vector and en¬ 
ergy, respectively. Eor simulations in a trapped cloud 
with N atoms, it is convenient to use the peak density 
of a non-interaeting Eermi gas in the Thomas-Eermi ap¬ 
proximation at the trap center, n' = {2jp)^jNmujjh^ 
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Figure 2: (color online). Upper panel: The phase jump 5(j) 
as a function of velocity in the non-topological phase. Lower 
panel: The corresponding mid-gap ABS energy levels. The ar¬ 
rows indicate the velocity at which the upper ABS (red cir¬ 
cles) touches the quasi-particle scattering continuum (shadow 
area), and the pair-breaking velocity Broyden’s method 
fails to find traveling soliton solutions when the energy of the 
lower ABS is close to zero. Parameters are as in Fig. [ 1 ] 


although the cloud itself is an interacting gas. We de¬ 
note the corresponding units with k'p and . Through¬ 
out this work, we consider only zero temperature. For 
trapped simulations we shall take the interaction param¬ 
eter 7 3, spin-orbit coupling strength Xk^/Ep = 1.5 

and an energy cut-off Ec = lOE^. Parameters for ho¬ 
mogeneous simulations are chosen to correspond to the 
relevant peak density of the interacting trapped gas. 

There are two different regimes for a ID spin-orbit- 
coupled Fermi superfluid [ 2 J, [23 , depending on whether 
the effective Zeeman field h is over a threshold he = 
-f (:^ Ep with our parameters for the trapped 
cloud). Once h > hc^ the superfluid becomes topolog¬ 
ically non-trivial and hosts Majorana solitons. Before 
presenting our main results on Majorana solitons, it is 
useful to understand how traveling solitons are affected 
by spin-orbit coupling in the non-topological phase. 

Non-topological phase. The spatial structure of the 
soliton order parameter in the nontopological phase {h < 
he) is illustrated in the lower panel of Fig. [T] for different 
soliton velocities. As the velocity increases, the dip in 
the order parameter profile becomes shallower, and its 
imaginary part develops structure and becomes larger at 
the soliton core. Consequently, the phase jump across the 
soliton decreases from tt, as shown explicitly in the upper 
panel of Fig. [2l This turn-to-grey procedure of traveling 
solitons has been predicted earlier both for BEC-BCS 
crossover superfluids 111 and BECs However, 

the presence of spin-orbit coupling leads to some inter¬ 
esting new features. 

The most striking feature is that the mid-gap energy 
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Figure 3: (color online). Time-dependent simulations of trav¬ 
eling solitons in a trapped, non-topological Fermi superfluid 
with h = 0.4E^ < he. The color represents the magnitude 
of density (in units of n'). By choosing the initial position 
xo, we generate two solitons, whose maximum velocity is 
O.lSu^ < Vh = 0 . 22 u^ (left panel) and 0.42u^ > Vh (right 
panel), respectively. The inset examines the universal relation 
for the soliton oscillation period with/without spin-orbit 
coupling at different interaction strengths (2.5 < 7 < 3.2). 
The period Ts from the time-dependent simulation is com¬ 
pared with TbcLG defined by the right hand side of Eq. ® 
and calculated from the time-independent BdG solutions. 


levels of soliton-induced ABSs now exhibit a pronounced 
velocity dependence, as seen from the lower panel of 
Eig. [21 Already at zero velocity, the ABS splits into two 
branches due to the combined effects of spin-orbit cou¬ 
pling and effective Zeeman field 26, 2^- With increasing 
the soliton velocity, the energy of the upper ABS grad¬ 
ually increases and merges into the quasi-particle scat¬ 
tering continuum at Vh — 0.22 vf, which is much smaller 
than the pair-breaking velocity O.dbui?. Any cou¬ 

pling between the upper ABS and the bulk continuum 
states [ 3 ^ then will destroy the soliton-induced ABSs 
and in turn make the soliton unstable. Thus, we an¬ 
ticipate that the soliton may decay when its velocity is 
beyond the threshold for example, by dissipating its 
energy in the form of sound waves. 

We have checked this conjecture by performing time- 
dependent simulations in harmonic traps, as reported in 
Eig. |3l By carefully selecting the position xq of the ini¬ 
tially stationary dark soliton, the maximum velocity Vm - 
reached when the traveling soliton passes the trap center 
- can be tuned. Eor Vm <Vh^ we find a stable oscillation 
of the traveling soliton (see the left panel). The oscil¬ 
lation period Tg seems to satisfy the elegant universal 
relation (see the inset). 


M* 

~w 


= 1 + 


hn d{S(l)) 
2M dva 


( 5 ) 


which was derived by treating soliton as a classical par¬ 
ticle liE|. Here, Tx = 27r /00 is the trapping period, M 
and M* are respectively the physical and inertial mass 
of the soliton and their difference is proportional to the 
derivative of the phase jump In contrast, at Vm > Vh^ 
the soliton gradually spreads out in the density profile 
and after a few periods we see only low-amplitude den- 
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Figure 4: (color online). Majorana soliton. The two upper 
panels report the time evolution of the magnitude \A{x,t)\ 
(left) and the phase (t){x,t) (right) of the order parameter for 
a Majorana soliton in a trapped topological Fermi superfiuid 
(h = 1.2E'p > he) with a maximum soliton velocity O.OTuf < 
Vh- The two lower panels show the density (left) and order 
parameter (right) of a Majorana soliton in the homogeneous 
configuration at different velocities with parameters 7 = 3.75, 
h = 1.71Ef and Xbr/Ep — 1.79. 


sity ripples (right panel). By examining the speed of 
these ripples, we identify them as sound waves. Our time- 
dependent simulations with spin-orbit coupling therefore 
indicate that the critical velocity of traveling solitons 
could be significantly smaller than Landau critical ve¬ 
locity Vp})^ which was found to be the relevant critical 
velocity without spin-orbit coupling liini. These result 
are still consistent, since without spin-orbit coupling Vh 
actually is close to the pair-breaking velocity [l2| . 

Topological phase. By increasing effective Zeeman field 
across he — Ep for a trapped Fermi cloud, the local 
energy gap (and hence the pair-breaking velocity) at 
the trap center closes and then re-opens. A topologi¬ 
cal superfluid emerges.The first sign of the existence of 
a velocity-independent Majorana soliton comes from the 
time-dependent simulations in harmonic traps, as shown 
in the upper panel of Fig. IH During the time evolution, 
the dip minimum in |A(x,t)| remains at zero and the 
phase jump 5(j){t) across the soliton is always pinned at 
TT (see also the inset in Fig. [5]). In the lower panel of 
Fig.m we check more rigorously the velocity dependence 
using Broyden’s approach. With increasing the soliton 
velocity in the topological phase, the density and pairing 
order parameter profiles remain essentially unchanged. 

To show the presence of Majorana fermions at the soli¬ 
ton core, we report in Fig. 5 the energy of the ABS as 
a function of the traveling velocity. Although in the co¬ 
moving frame the energy E^^g increases (linearly) with 
the velocity, the energy in the laboratory frame, 
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Figure 5: (color online). The ABS energy of the Majorana 
soliton as a function of the soliton velocity in the co-moving 
frame (upper panel) or in the laboratory frame (lower panel). 
We note that in the topological phase, the number of the ABS 
states decreases to one [^, if we count only positive energy 
levels. The inset examines the 7 r-phase of Majorana solitons. 
Parameters as in Fig. [H 


which is related to the co-moving energy by 


77'Lab 77' 

AUQ = ^ 


_ iMov 

ABS — ^ABS 




I (6) 


is precisely zero [^. This is expected behavior for a 
Majorana fermion, which must have zero energy due to 
the particle-antiparticle symmetry. Together with the 
observed continuity with the zero velocity case 0113, 
we conclude that the moving soliton in the topological 
phase indeed hosts Majorana fermions.The properties of 
the Majorana soliton at finite velocity can be made plau¬ 
sible from the universal relation if we assume its 
validity in the topological phase. We recall that the den¬ 
sity notch in Majorana solitons i^absent [ 2 ^, 271 and 


hence the physical mass vanishes Equation (j^ im¬ 
mediately implies that the derivative of the phase jump 
is zero, since the oscillation period should be finite. This 
leads to a constant tt phase jump, irrespective of the soli¬ 
ton velocity. In turn, the magnitude of order parameter 
should vanish at the soliton core. 

It is worth noting that Eq. m can hardly be used to 
predict the oscillation period of Majorana solitons, as the 
ratio between zero mass and zero derivative of the phase 
jump is undetermined. How to amend the universal rela¬ 
tion for Majorana solitons needs further exploration. Eor 
a large velocity, we find a similar situation as in the non- 
topological case (see the upper panel of Eig. [5j). Naively, 
we anticipate that Majorana solitons may cease to ex¬ 
ist once Vs > Vh- However, time-dependent simulations 
in traps suggest that the instability via emitting sound 
waves occurs at a very long time scale, presumably due 
to the weak coupling between the ABS and quasi-particle 
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continuum. This again indicates the robustness of Majo- 
rana solitons. 

Experimental observation of Majorana solitons. The 
proposed experimental scheme in Fig. [1] is easy to set up 
[ 22 I, l28|, although the realization of ID topological su¬ 
perfluids is difficult due to the lack of efficient cooling 
techniques [22|. Majorana solitons can be created by im¬ 
printing a sharp phase jump [mill. The observation of 
their oscillations seems to be an experimental challenge, 
as the density profile of Majorana solitons remains flat 
One then has to measure the local pairing order 


parameter. A suitable detection technique is the spa¬ 
tially resolved radio-frequency spectroscopy [s^, which 
may give information about the local order parameter, 
provided that the size of Majorana solitons is compara¬ 
ble with the spatial resolution of spectroscopy. 

Conelusions and outlook. We have predicted the ex¬ 
istence of an exotic member to the soliton family - the 
Majorana soliton - which may exist universally in any 
ID topological superfluids including p-wave superfluids 
and semiconductor/superconductor nanowire structures 
[ 35 I I. Our results may also be applicable to 2D topolog¬ 
ical superfluids. In that case, it would be interesting to 
examine the possibility of finding a moving vortex that 
is able to host a single Majorana fermion in the vortex 
core. In analogy to Majorana solitons, the properties of 
such a vortex would be insensitive to its velocity. 
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